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We estimate generic statistical properties of a structural credit risk model by considering an ensemble of 
correlation matrices. This ensemble is set up by Random Matrix Theory. We demonstrate analytically that 
the presence of correlations severely limits the effect of diversification in a credit portfolio if the correlations 
are not identically zero. The existence of correlations alters the tails of the loss distribution considerably, 
even if their average is zero. Under the assumption of randomly fluctuating correlations, a lower bound for 
the estimation of the loss distribution is provided. 



I. INTRODUCTION 



The financial crisis of 2008-2009 clearly revealed that an improper estimation of credit risk can lead to dramatic 
effects on the world's economy. The vast underestimation of risks embedded in credits for the subprime housing 
markets induced a chain reaction that propagated into the worldwide economy. A better estimation of credit risk 



(see, eg,|Bluhm et al.| ( |2002[ ); [Bielecki and Rutkowski | ( |2005[ ); [Duffie and Singleton] p003| ); |Lando| p004| ); [McNeil 
et al. ( 2005p ) is therefore of vital interest. 



We can distinguish two fundamentally different approaches to credit risk modeling (see, eg, Giesecke (2004|): the 



and Scholes ( 


19731 and 


Merton 


(1974 



to maturity. The equity of the company is modeled by a stochastic process, for example describing the price of the 
company's stock. Thus, it can be seen as a creditor's call option on the obligor's assets. The risk of default and the 
associated recovery rate, the residual payment in case of a loss, are modeled by the company's equity at maturity. 
Reduced-form models attempt to capture the dependence of default and recovery rates on macroeconomic risk factors. 
Both quantities are modeled as independent stochastic variables. Some well known reduced-form model approaches 



can be found in Jarrow and Turnbull (1995); Jarrow et al. (1997); Duffie and Singleton (1999); Hull and White 



(2000); Schonbucher (2003). First passage models were first introduced by Black and Cox (1976). Although they 
are often referred to as structural models, they do, in fact, follow a mixed approach. Similar to Merton's model, 
the market value of a company is modeled by a stochastic process. However, in the first passage models a default 
occurs whenever this market value hits a certain threshold for the first time. The recovery rates are typically modeled 
independently, for example, by a reduced-form model, see Asvanunt and Staal ( 2009a|b" ) , or are even assumed to be 
constant, see, eg, Giesecke (2004). Recent approaches aim at improving fist passage models by including the chance of 
full recovery, even if a company's market value is below the threshold, see Katz and Shokhirev (2010), and estimating 
correlations between default probabilities of industry sectors, see Rosenow and Weissbach (20091. Reduced-form and 
first passage models are implemented in commercial software solutions, for example, CreditMetrics by JP Morgan 
(see Gupton et al. (1997)), CreditPortf olio View by McKinsey & Company ( |1998p or CreditRisk+ by Credit Suisse 
(1997). As there can be a strong connection between default risks and recovery rates, the chances of large losses are 



often underestimated in the reduced-form and first passage models, see Schafer and Koivusalo (2011 ); Koivusalo and 



Schafer ( 2011, ). Structural models do not require this separation. 



Structural models provide a "microscopical" tool to study credit risk as the defaults and recoveries are traced back 
to stochastic processes modeling the state of individual obligors. For a portfolio of credits, such as collateralized debt 
obligations (CDOs), correlations represent a key factor that influences its risk. The benefit of diversification, ie, the 
reduction of risk by increasing the portfolio size, is severely limited by the presence of even weak correlations. This 
has been demonstrated for the case of constant positive correlations, both in the first passage model with constant 
recovery (Schonbucher ( [2001 ); Glasserman (2003)) and in the Merton model (Schafer et al. (20071; Schafer and 



Koivusalo (2011 )). The key problem in estimating the credit risk of a realistic portfolio is of course the huge number 



of parameters involved. This is precisely where approaches from statistical physics can be most helpful: the state of 
a system with many degrees of freedom is, under certain conditions, described by few macroscopic observables. In 
the thermodynamic equilibrium, these are energy, temperature, pressure, etc. Ergodicity holds, ie, time and ensemble 
average yield the same results. A somewhat similar situation exists for spectral statistics in quantum chaotic systems. 
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Variable Description Unit 

K Number of assets - 

T Time of maturity [year] 

CTfc volatility of the fc-th asset [year]"^''^ 

Hk drift of the fc-th asset [year] ~ ^ 

N Parameter to control correlations, - 

N —> <x: uncorrelated limit 

Vkfi start price of the fc-th asset [currency] 

Fk face value of the fc-th asset [currency] 

TABLE I: Input of the structural credit risk model. 



see Guhr et al. (1998). A moving average over one long spectrum equals an ensemble average over random matrices, 
if the number of levels is very large. Originally, random matrix theory was developed in the 1950s to describe the 
spectra of heavy nuclei, see Mehta (2004 1. Here we transfer this idea to large credit portfolios involving correlated 
assets. In the case of a great many contracts, we expect a self-averaging property which then should allow to average 
over an ensemble of random correlation matrices. We manage to carry out this approach largely analytically. We 
obtain estimates for the price and the loss distribution in which the complicated effects of all correlations are indeed 
reduced to a single parameter measuring the correlation strength. This paper is organized as follows. A structural 
credit risk model with random correlations is considered in section |Il] We demonstrate an application of this model 
in section 



III and conclude our findings in section IV 



II. A STRUCTURAL CREDIT RISK MODEL 



Our model is based on Merton's original model. Our aim is to analytically describe the impact of correlations on 
the losses of a credit portfolio. A default occurs if the price Vk of the fc-th asset is below the face value at maturity 
time T. The size of the loss then depends on how far Vk is below the face value Fk- We assume that the prices in a 
portfolio of K assets follow a geometric Brownian motion. An overview of the model's input parameters is given in 
Tab. H 



A. Average price distribution 

For the sake of simplicity, let us first consider the case of a Brownian motion for the average price distribution. 
Later on this can be easily mapped to the geometric Brownian motion by a simple substitution. For a Brownian 
motion, the probability density function (pdf ) of the price vector V of K assets at maturity T is described by 

P<»-.(v,E) . -^^=p(-±(r-.r,.E-.(v-,,r)) ,1) 

Here, E is the covariance matrix and fi is the drift vector. For later convenience we can express this as a Fourier 
transform. 



S) = J exp {-tLu\V - ^iT)) exp 



-!^a;^Sw ) d[a;] (2) 



Equation ([T]) gives the pdf of asset prices in the case of a correlated Brownian motion. However, we are not interested 
in the impact of a specific correlation matrix. Instead we want to estimate the general impact of correlations. To 
this end, we want to average over all possible correlation matrices and disclose the general statistical behavior of the 
system. This will enable us to make a profound statistical statement. 

We use a random matrix approach to calculate the average price distribution for random correlations where the 
average correlation level is zero. By averaging over all possible combinations of random variables, we obtain the 
average price distribution {p^'^^\V)) under these assumptions. To achieve this we replace the covariance matrix S by 



= SWW^'S 



(3) 
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where S = diag(cri, . . . , aK) contains the standard deviations and W € M^^^ is a random matrix. The entries of W 
are independent and Gaussian distributed, 



(corr) 



KN 



exp 1 -ytr 



(4) 



with variance This results in a correlation matrix WW^ with average correlation zero. With the parameter N 

we can control how strongly the entries of WW^ fluctuate. For — > oo, we obtain the pdf of a unit matrix. This 
represents the uncorrelated case. For N > K, we obtain an invertible covariance matrix with random entries. The 
case N < K is disregarded as the resulting matrix is not invertible which is usually required for applications in risk 
management. When inserting this ansatz into Eq. ([2]), we obtain 



= / p'^''°"\w)p^"'''\V, SWW'< S)d[W] 

-NK 



l2 



^ f ( 

^ J (- 



1 



^/det{NI + TSuju^S) 



■N 



(5) 
(6) 



where / denotes the unit matrix. A detailed derivation is given in appendix [X] Here we choose ^ — 0. We will 
reintroduce the drift later on, when we make the substitution for the geometric Brownian motion. The determinant 
can be written as 



det(iV/ + TSujuj^S) = (l + ^uj'fSSiu) 
because the matrix Sujlu^ S has rank one. Hence, we arrive at 

1 



This integral can be calculated by using the Gamma function (see Olver ( 1974[ )) in the form 



(7) 



(8) 



Fix) 



exp {—az) dz , x > 0, a > 



(9) 



We identify a"^ with ((1 + {T/N)uj^ SSlu))-^/^ and obtain 



{2n)i< T{N/2) 

oo 

X / z(^^^^ exp (— z) 




K 



exp 



' 4Tz ^ 0-2 

fc=i , 



dz 



(10) 



as worked out in appendi x [B| T his integral is a representation of the Bessel function of the second kind /C of the order 
{K - N)/2, see |Watson| ( |1944| ). Thus, we obtain 



(p(--)(y)) 



N 2^- 



2t:T T{N/2) 



\ 



T 



fe=i 



(11) 



for the average distribution of {V) if assuming a randomly distributed correlation matrix and an average cor- 
relation level of zero. We stated earlier that we include N in the distribution of the random matrices W in order 
to render the variance of the average price distribution A'^- independent. The variances only depend on T and Cfe, as 
discussed in appendix [C] The parameter iV is only used to control the correlations. In hyperspherical coordinates. 
Equation 11 depends only on the hyperradius 



(12) 
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This leads to the expression 



N 2i- 



2ttT T{N/2)' 



N 
T 



(13) 



for the hyperradial density function, cf. appendix [C| We illustrate this density function in Figure [T] for K — and 
different values of N . 

We obtain the average price distribution in case of a geometric Brownian motion by a simple substitution — > V^, 



1 i 



1-kT 



^ 1 



r(iV/2) ^11.,^,;^ T^^ a 



k 



T 



E 



-1 '^fe 



fc=i 



with 



14 = In 



14.0 



(14) 



(15) 



Here, the parameter refers to the standard deviation of the underlying Brownian motion, ie, the volatility of asset 
returns. The resulting prices thus have the variance 



(Tfc = ^exp {alT + 2m) (exp (a^T) - \)Vl^ (16) 

where Vfc^o are the starting prices at i = 0. Figure [2] shows the distribution of prices based on a geometric Brownian 
motion, as given in Eq. ( 14 ) . The findings are similar to the case of the Brownian motion. While we obtain a narrow 
but heavy-tailed distribution for N = K, the distribution slowly approaches an uncorrelated bivariate log-normal 
distribution with increasing values of N . 



B. Loss distribution 



We now turn to the calculation of the loss distribution. A default occurs if the price 14 at maturity T, corresponding 
to the obligor's equity, is below the face value Fk- The size of the loss is given by the difference of Fk and Vk- Even 
if a loss occurs, the creditor might not lose all money that he lent, because the obligor is still able to pay back the 
amount Vk- In order to compare losses in a portfolio of credits, we have to normalize them by the corresponding face 
value. We define the normalized loss Lk of the fc-th asset as 

^ , Vk<Fu (defauh) ^^^^ 
, else (no default) 
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FIG. 2: Illustration of the average price distribution {p^™^\V)) with a geometric Brownian motion for T — 1, K — 2, 
Vfc.o — 100, ^ — 0.05 and different values for N. Both distributions have the identical standard deviation ct « 16 
(cr = 0.15). For = 2, we obtain a heavy-tailed distribution while the uncorrected limit is reached for N = 100. 



We observe that the prices have to be positive in Eq. (17). Therefore we assume in all further considerations that the 



underlying process of the price distribution follows a geometric Brownian motion. 

When calculating the overall loss of a portfolio, we have to weight each loss by its face value in relation to the sum 
of all portfolio face values, 



K 



Fk 



k=l 



(18) 



We integrate over the pdf of prices and filter for those that lead to a given total loss L. By the above stated definitions, 
we can define a filter for the total loss at maturity time T. In the next step we express the filter using a Fourier 
transformation. Eventually, we separate those terms that correspond to a default and those that describe the price 
above the face value Fk. 



(loss) 



exp —ivL + iv 



Jd[V]i"^-\V)5 (^L~f2fkLk^ 

f^hLk] 

k=i / 

^- I dv exp {-ivL) I d[V] exp I iv fkLk ] P*'""' (V) 

-L V fc=i / 

1 r 

— / dv exp(—ivL) 
In J 



(19) 
(20) 
(21) 



K 

n 

k=l 



dVk exp iiyfk 1 - 



Vk 
Fk 



dVk 



(22) 



Here, the expression in the square brackets acts as an operator, because p^™^\V) does not necessarily factorize. We 
will use this ansatz to calculate the average loss distribution in the next section. However, Eq. (22 1 can be used 



to calculate the loss distribution if the actual price distribution is known, ie, the statistical dependence and the 
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underlying process are estimated. To prepare for this, it is handy to write Eq. (22) as a combinatorial sum 



(loss) 



+00 

— / dveyi'p{—wL) 
2vr J 



(23) 



xEE 

k=l3 = l 



n 



IGPcrm(j,A;,K) 



dV; exp ivfi 



Fi-V, 



F, 



n 



qe{i...K} 

\Porm(3,fc,K) 



n(mv) 



iV) 



where Perm(j, fc, K) is the j-th permutation of k elements of the set {1 . . . K}. For example, ii K — 3 and k — 2, 
we obtain, Perm(l,2,3) = {1,2}, Perm(2,2,3) = {2,3} and Perm(3,2,3) = {1,3}. However, Eq. (|24l ) mig ht need 



to be estimated numerically, depending on the complexity of the price distribution p(™^-'(V^). In section II D we will 



simplify this combinatorial sum for a homogeneous portfolio and the average price distribution {yTT- 



C. Average loss distribution 



Now we have developed all necessary tools to model the average distribution of losses, under the assumption of 
random correlations and an average correlation level of zero. We start by inserting the average price distribution in 
a component- wise notation (see appendix ^ into the loss distribution ( 22 ) , 



(loss) 



1 



+00 



27rr(A^/2) 



exp (— z) / dv exp {—ivL) r{i', z) 



with 



K 



dVfe exp ii^fi 



k=l 



X exp 

2akVkVTTzT 



Fk-Vk 
' Fk 



N{\n{Vk/Vu^o)^{^i-a^/2)Tf 



4:zT(jI 



We carry out a second order approximation of this expression in appendix [D] and arrive at 



(loss) 



1 



27rr(iV/2) 



dz z( 2 1) exp (— z) 



1 



exp 



(L-TOl(z))^ 

2m2{z) 



with 



K 



mi(z) = ^ fkmi,k{z) 
fc=i 

K 

m2{z) = ^ /|(m2,fc(z) - mi^k{zf) 



fe=i 



and 



N 



2akVTr^ J Vk 



1 fFk-VkV 



Fk 



X exp 



N{\n{Vu/Vk,o) -{fi- a^/2)Tf 
AzTal 



dVk 



(24) 



(25) 
(26) 

(27) 
(28) 



(29) 



However, the convergence radius of the power series expansion involved in this approximation is one. Although we 
consider large portfolios ie, is small, v runs from — cxi to +00. This second-order approximation might describe 
the default terms adequately. However, the non-default terms, corresponding to a delta peak at L = require v to 
run from —00 to -l-oo. Thus, the non-default terms cannot be approximated using this second-order approximation. 
To circumvent this problem we develop an improved approximation in section [ll D| 

Due to the complexity of rhi{z) and TO2(z), the z integral needs to be evaluated numerically. We present this for 
the example of a homogeneous portfolio. 
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D. Homogeneous portfolios 

In case of a homogeneous portfolio, in which all credits have the same face value F}, = F and the same variance 
af. = and initial value Vfc_o = Vb, the weights can be simplified to 

h = ^ (30) 

As mi^kiz) and mi^kiz) become identical for every k, we denote them by toi(z) and toi(z) leading to 

mi(z) = mi{z) (31) 
m2(z) = -^(m2(z) - mi(z)^) (32) 



mj{z) 



F 

N f I fF-V^^ 



2aVT^ J y \ F 





^2 /o^T^^2 



x..p|- ^"""^/^°'j7'''/^>^' |dV (33) 



Here y is a scalar and we only have to calculate a single integral over V. After inserting this into Eq. (26), we can 
calculate the loss distribution for a homogeneous portfolio in the second order approximation. 



E. Improved approximation for a homogeneous portfolio 

The second order approach can be improved by approximating the individual terms of the loss distribution instead 



of approximating the expression as a whole, similar as discussed in Schafer et al. (2007). In case of a homogeneous 



portfolio the combinatorial sum in Eq. (24) reduces to 

+ C30 



J=0 

with the non-default term {r'^^^)\^ where 



P^'°'''\L)) = ^^^j;^ / dz z(^-i)exp(-z) j dve^v{-iuL) 

-oo 

X:('^)(r(°'(-,-))'(^'^°'(^))''" (34) 



r(^°) = / d^^ 



^aVyJuzT 

F 

2 /o^T^^2 



1^1^ / V^(ln(FA«-|.-,V2)T) 



and the default term (r^^^iy, z))"* where 

F 

A^\v,z)^ f dFexp^*^^"^ 



K F 

N ^ ( N{HV|V,)-{^,-a^/2)Tf \ 

71^2 (37) 
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In the homogeneous case V is a scalar variable. The approximation follows the same principles as in the previous 
section, resulting in 



+00 



diy exjp {—ivL) (^A^^v, z) 



+00 



dj/exp ( iiy ( — mi(z) — L 
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{m2iz) - mi(z)'' 



j (m2(z) - mi(z)2) 



exp 



2^2 

iLK-jmi{z))^ 

2j (to2(z) - 



(38) 
(39) 



In this approximation, the non-default terms given by Eq. ( 36 ) can be calculated exactly. They correspond to a delta 



peak at L = 0. Another advantage over the approach presented in Eq. (261 is that the approximation is performed for 



each number of defaults j separately and weighted by j /K accordingly. Here, the omitted third term is of the order 
j / and thereby much smaller than the third term of the simple second order approximation (33), which would be 
of the order Thus, when approximating each term in the combinatorial sum separately, we obtain an improved 

result. Insertion into (34) leads to 



27rr 



^i:(-)/a=.(»-0»p,-, 

{LK~]m^{z)f 



27rX2 



j {m2{z) - mi{zy 



exp 



2j (to2(z) - TOi(z)2) 



HF/Vo) - (Ai - ay2)T 



2aVzT 



(40) 



which is the final result. 



III. APPLICATION 



We now apply the analytically developed model to a specific example. To analyze the impact of correlations, we 
calculate the loss distribution for different homogeneous portfolios with sizes K = 10, K — and K = 100 with the 
parameters Vb = 100, /i = 0.05, a = 0.15, / = 75 and T = 1. As stated in the previous section, we can control 
the amount of correlation in our model with the parameter N . Since we only consider correlation matrices with full 
rank, we obtain the strongest correlations if we choose N = K. For N ^ 00, the correlation matrix becomes the 
unit matrix. Thus, this represents the transition to a system without correlations. As we have to evaluate the loss 
distributions numerically, the limit iV — >■ 00 has to be properly interpreted. We need to identify a value for which this 
convergence is valid in good approximation. Figure [3] illustrates the loss distribution for if = 10 and different values 
of N. Our study indicates that a value of = 30Ar is a good choice for approximating the uncorrelated case and is 
still numerically feasible. The results are presented in Fig. |4j For all portfolio sizes, K = 10, K = 50 and K = 100, 
we obtain heavier tails of the loss distribution of the correlated portfolio compared to the uncorrelated case. Even the 
simple approximation, represented by the dashed blue curve, exhibits these heavy tails. With the inserted logarithmic 
plots, we can identify a nearly power-law decay of the loss distribution for the correlated case. 

The distributions become narrower for larger values of K. However, the tails of the correlated case remain heavier 
than those of the uncorrelated case. While both approximations yield similar results for K = 10, their difference 
becomes larger with K. As both approximations have to be performed numerically, the improved approximation is 
always favored. However, the tail behavior remains the same, even for the simple approximation, as indicated by 
the logarithmic scaled inserts in Fig. |4j This is a strong indication that the tails of the loss distribution are vastly 
underestimated if correlations are not taken into account. 

Due to the approximation, the normalization of the loss distribution is not exact. Especially the normalization of 
the simple approximation is problematic for large values of K. The normalization might also be used as an indication 
for the quality of the approximation. The improved approximation exhibits a delta peak at L = 0, as the non-default 
terms can be calculated exactly. However, the interval [0;0.0002[ was not evaluated due to numerical feasibility. 
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0.000 0.002 0.004 0.006 0.008 0.010 



L 

FIG. 3: The loss distribution for K ^ 10, a = 0.15, /i = 0.05, T = 1, Vo = 100, / = 75 and different amounts of 
randomness in the correlation matrix. N ^ K (solid black), N = 2K (dashed blue), N — lOK (dotted red), 

N = 30K (dot-dashed green). 



In our example, we do not vary the maturity time T, ie, we choose T — 1. One can increase T to estimate the 
evolution of the loss distribution. However, this evolution depends strongly on the drifts fik and standard deviations 
CTfc. Depending on their value, the exposure to default risk can either increase or decrease. 



IV. CONCLUSIONS 



To assess the risk of a credit portfolio, it is crucial to take correlations between obligors into account. We consider 
the Merton model, in which defaults and recoveries are determined by the underlying asset processes. The correlation 
matrix of the asset returns has to be estimated from historical time series. This is not always easy, because the 
correlations change in time, ie, they are non-stationary. Since only time series up to a certain length can be used, 
the correlation coefficients contain a specific type of randomness, see Laloux et al. (19991. Several methods have been 
put forward to estimate and to reduce this "noise". Thus, we assume that such a noise reduction has been done. The 
corresponding "true" correlation coefficients and matrices are the proper input for the structural credit risk model of 
the Merton type that we consider. We discussed this issue of noise reduction to emphasize that the random matrix 
approach in that context is based on a very different motivation as compared to the present application. 

Searching for generic properties, we devised the present random matrix approach. Instead of calculating the portfolio 
loss distribution for a specific correlation matrix, we average over an ensemble of random correlation matrices. Our 
approach transfers concepts of statistical physics. In quantum chaos, the average over an individual, long spectrum 
equals the average over an ensemble of random matrices, if the level number is very high. We expect that a similar 
self-averaging property also holds here. This line of reasoning is supported by the following consideration: The 
correlation coefficients are varying functions in time, because the business relations of the companies change. This 
implies that a correlation matrix over a somewhat longer period in time is a varying quantity, ie, it corresponds to 
some kind of ensemble. 

In our model the average correlation level is zero and we assume that there is no branch structure in the correlations. 
The fluctuation strength of individual correlations is controlled by a single parameter. This ansatz allowed us to 
estimate generic statistical properties of the Merton model. Some features are not taken into account which are 
present in empirical data, such as jumps or an overall positive correlation level. Those features are difficult to treat 
completely analytically. However, even in our simple setup we obtain a heavy-tailed loss distribution. In this sense 
our model can be used to estimate a lower bound of the risk embedded in a credit portfolio. 

Our results clearly demonstrate that the risk in a credit portfolio is heavily underestimated if correlations are 
not taken into account. Even for random correlations with an average correlation level of zero, we observe very 
slowly decaying portfolio loss distributions. In contrast, the probability of large losses in uncorrelated portfolios is 
significantly reduced within the Merton model. 
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(a) K=10 
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(b) K=50 
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(c) K=100 

FIG. 4: The loss distribution of a homogeneous portfoho with a — 0.15, /i — 0.05, T — 1,Vq — 100, f — 75 and 
different values of K. The blue dashed line represents the simple approximation; the solid black line represents the 
improved approximation. Both have been calculated for the strongest random correlations, N = K. The 
uncorrelated case is given by the red dotted line, calculated with the improved approximation with N = SOK. 
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The results are especially relevant for CDOs, bundles of credits that are traded on equity markets. CDOs are 
constructed in order to lower the overall risk. The components of a CDO can be exposed to large risks. It is often 
believed that the CDO has a significantly lower risk. We showed that this diversification only works well if the 
correlations in the credit portfolio are identical to zero. 
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Appendix A: Identifying a Gaussian integral in average the price distribution 



(p(«i-)(y))= I p<^''°"\w)p''"''-\V, SWW^ S)d[W] (Al) 

^ \ 2K+KN 



IV^^y" exp [-^irWW'^^ 
X j exp {-iuj^V) exp (^—uj'^ SWW^ Su^ d[uj]d[W] (A2) 

y-^j J exp {-iuj^V) 

X J exp J^-^trWW^I'^ 

X exp ( - ^^wl'SWW^t^j^^ d[Ty]d[w] (A3) 



In the last steps, we took advantage of the fact that the term co^ SWW^ Sco is a scalar, which can evidently be written 
as trace. As the trace is invariant under cyclic permutation, we can express this term as tT{WW^ Scoco^ S) . Hence, we 
write 



/ 1 \2K+KN I- 

X J cxp(^~^tr{WW\NI + TSojio^S))jd[W]d[io] (A4) 

/ , X 2K+KN . 

= [-^) ^ /expH.t^) 

X y exp 1^-^ ^ {Wl{NI + TSojJS)wS^ dK]d[w] (A5) 

J exp(-^w''{NI + TSojJS)w)^d[w]] d[uj] (A6) 

where / is the unit matrix. The last step can be accomplished, as the components of W are independent identically 
distributed, hence we can denote the n-th cohinin vector of W, by w. Thus, we can simpHfy the integration over 
the matrix W to the integration over the vector w G to the power of A''. The integral over d[w] is simply a 
Gaussian integral, as indicated by Xj, Xj and Aij. As Wn consists of K components, this gives an additional factor 

ZTT and thus leads to 
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Appendix B: Derivation of component- wise average price distribution 
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Appendix C: Normalization of the average price distribution 



We consider the variance of the z-th price Vi^ 

var(F,) = y" d\V\ V^lJ--^^\V)) 
We can solve this integral by using hyperspherical coordinates. 
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for a chosen component Vi/ai. Now we can write the integral in Eq. (CI I as 
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15 



with the surface of the corresponding i^-dimensional sphere 



We obtain 



2 T T{N/2 + 1) 2^ 



Thus, the variance of every Vi only depends on the standard deviation Ui and the time T. 



Appendix D: second order approximation of average loss distribution 



We rearranged the constants so that each term in r{v, z) is normaUzed to unity. The quantity r{v, z) can now be 
written as 

K r / / ^ ^ OO \ 



r{v^ z) = exp j In 
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dVfe exp ivfk 1 
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We expand q (i/, Fk) as the power series 



(Dl) 



(D2) 



Due to the nor mahzation of (p^'"^^^)): after insertion into Eq. ( |Dl[ ) the non-default term and the integral over first 
term of Eq. (D2 1 become one. Thus, we can start the sum at j = 1 and obtain 



r(i., z) = exp (^t^ln (^1 + g iii^m,- ^(z)^ ^ 



(D3) 



with 
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(D4) 



The integrals in Eq. (D4) can be expressed with the generalized hypergeometrical function pFq. However, the integral 
representation (D4) is more intuitive. Moreover, for explicit m = 1, 2 the integrals can be calculated in a closed form, 
although this results in bulky expressions. 

Expanding the logarithm and collecting all terms up to the second order in fk yields 
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Now we can solve the v integral leading to Eq. (26). 



